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Abstract 

We study the matching of the next-to-leading logarithmic approximation (NLLA) onto 
the fixed next-to-next-to-leading order (NNLO) calculation for event shape distributions in 
electron-positron annihilation. The resulting theoretical predictions combine all precision 
QCD knowledge on the distributions, and are theoretically reliable over an extended kinc- 
matical range. Compared to previously available matched NLLA+NLO and fixed order 
NNLO results, we observe that the effects of the combined NLLA+NNLO are small in the 
three-jet region, relevant for precision physics. 
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Event shape distributions in e + e~ annihilation processes are classical hadronic observables 
which can be measured very accurately and provide an ideal proving ground for testing our 
understanding of strong interactions. The deviation from simple two-jet configurations, which 
are a limiting case in event shapes, is proportional to the strong coupling constant, so that 
by comparing the measured event shape distribution with the theoretical predictions, one can 
determine the strong coupling constant a s . At LEP, a standard set of event shapes was studied 
in great detail [1—4]: thrust T [5] (which is substituted here by r = 1 — T), heavy jet mass p [6], 
wide and total jet broadening B\y and Bt [7], C-parameter [8] and two-to-three-jet transition 
parameter in the Durham algorithm Y3 [9]. The definitions of these variables, which we denote 
collectively as y in the following, are summarised in [10]. The two-jet limit of each variable is 
y^O. 

The theoretical prediction is made within perturbative QCD, expanded to a finite order in 
the coupling constant. This fixed order expansion is reliable only if the event shape variable is 
sufficiently far away from its two-jet limit. In the approach to this limit, event shapes display 
large infrared logarithms at all orders in perturbation theory, such that the expansion in the 
strong coupling constant fails to converge. Resummation of these logarithms yields a description 
appropriate to the two-jet limit. To explain event shape distributions over their full kinematical 
range, both descriptions need to be matched onto each other. Until very recently, the theo- 
retical state-of-the-art description of event shape distributions was based on the matching of 
the next-to- leading- logarithmic approximation (NLLA, [11]) onto the fixed next-to- leading or- 
der (NLO, [12-14]) calculation. Using the newly available results of the next-to-next-to-leading 
order (NNLO) corrections for the standard set of event shapes [15-18] introduced above, we 
derive here matching of the resummed NLLA onto the fixed order NNLO. 

For two-particle final states, all above event shape variables have the fixed value y = 0, con- 
sequently their distributions receive their first non-trivial contribution from three-particle final 
states, which, at order a s , correspond to three-parton final states. Therefore, both theoretically 
and experimentally, these distributions are closely related to three-jet production. 

Fixed-order QCD corrections to event shape distributions were calculated long ago to next- 
to-leading order (NLO, [12-14]), and most recently to next-to-next-to-leading order (NNLO, [15- 
17]). At a centre-of-mass energy Q and for renormalisation scale ft, they take the form: 

Zr(v,Q,n) = as(^)^(y) +al{fi)^(y,x^) +a 3 s (fi)^-(y,x^ + 0(af) , (1) 



Chad dy dy s dy s dy 

where 

a x -t (2) 

as ~ 2vr ' Q ' ( } 

and where A, B and C are the perturbatively calculated coefficients at LO, NLO, NNLO, 
normalised to Ohad; explicit relations are given in [17] (note the different convention for the Pi 
coefficients). 

The resummation of large logarithmic corrections in the y — > limit starts from the inte- 
grated cross section: 

Sfeft^-Lf^,, (3) 

Chad JO dX 
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which has the following fixed-order expansion: 

R{y,Q,l*) = 1 + A (y) a s (p) + B (y, x^) a 2 s (fi) + C (y, x M ) a 3 s (p) . (4) 

The fixed-order coefficients A, B, C can be obtained by integrating the distribution ([1]) and using 
-R(ymax, Q, fJ-) = 1 to all orders, where y ma x is the maximal kinematically allowed value for the 
shape variable y. 

In the limit y — > one observes that the perturbative a"-contribution to R{y) diverges like 
a™L 2n , with L = —\ny (L = — ln(y/6) for y = C). This leading logarithmic (LL) behaviour is 
due to multiple soft gluon emission at higher orders, and the LL coefficients exponentiate, such 
that 

In R(y) ~ Lgi{a s L) , 

where g\(a s L) is a power series in its argument. 

For the event shapes considered here, leading and next-to-leading logarithmic (NLL) correc- 
tions can be resummed to all orders in the coupling constant, such that 

R (y, Q, //) = (1 + da a ) e (^iK^KiO) , (5) 

where terms beyond NLL have been consistently omitted, and \x = Q (xn = 1) is used. In the 
case of the C-parameter further large logarithms around C m 0.75, produce a so-called Sudakov 
shoulder in the distribution due to soft gluon divergences within the physical region [24]. 

By differentiating expression ([5|) with respect to y, one recovers the resummed differential 
event shape distributions, which yield an accurate description for y — > 0. The first complete cal- 
culation of next-to-next-to-leading logarithmic (NNLL) corrections to event shape distributions 
is available for the energy-energy correlation function [27], which is not part of the standard 
set of event shape observables. The application of soft-collinear effective field theory to event 
shape distributions [28] promises to yield results beyond NLL. Most recently, this formalism was 
applied to compute the resummed thrust distribution beyond NLL accuracy [29]. 

Closed analytic forms for the LL and NLL resummation functions g\(a s L), g2{a s L) are 
available for r [19], p [20], B w and B T [21,22], C [23] and Y 3 [25]. For the convenience of the 
reader, we collect them in uniform notation in an Appendix. They can be expanded as power 
series, such that: 

oo i+1 

\nR(y, Q,fi) = J2J2 Gi^-n^V^ , (6) 

i=l n=l 

To obtain a reliable description of the event shape distributions over a wide range in y, it 
is mandatory to combine fixed order and resummed predictions. To avoid the double counting 
of terms common to both, the two predictions have to be matched onto each other. A number 
of different matching procedures have been proposed in the literature, see for example [10] for 
a review. The by-now standard procedure is the so-called In i?-matching [11]. In this particular 
scheme, all matching coefficients can be extracted analytically from the resummed calculation, 
while most other schemes require the numerical extraction of some of the matching coefficients 
from the distributions at fixed order. Since the fixed order calculations face numerical instabil- 
ities in the region y — > 0, these matching coefficients can often be determined only within large 
errors. We shall therefore consider only the In i?-matching here. The In ii-matching at NLO is 
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described in detail in [11], where the authors also anticipated the fixed-order NNLO corrections 
to be available shortly, and briefly outlined this matching scheme to NNLO. 
In the In /^-matching scheme, the NLLA+NNLO expression is 

In (i? (y, as)) = Lg x (a s L) + g 2 (a s L) 

+ a s U(y)-G 11 L-G 12 L 2 



+ a| (b (y) - \a 2 (y) - G 22 L 2 - G 23 L 3 

+ 4 (c (y) -A{y)B (y) + ±A 3 (y) - G 33 L 3 - G U L^ . (7) 

The matching coefficients appearing in this expression can be obtained from ([6]) and are listed 
in Table [TJ In the matching of Y 3 , the constants Ti depend on the jet algorithm [25], in general, 
they can be determined only numerically. For the Durham-algorithm, one finds J- 2 = — 7r 2 /32 
and T 3 = 0.0620 ± 0.0100 [30], using the semi-numerical resummation method described in [26]. 
Numerical values of the matching coefficients for N = 3, Np = 5 are given in Table [2J 

To ensure the vanishing of the matched expression at the kinematical boundary y max , the 
further substitution [10] is made: 

L _ £ = i ta ((j!Ly_(_!5_y +1 '), (8) 

P WXLVJ XXLVm^J J 

where yo = 6 for y = C and yo = 1 otherwise, p = 1 and xl = 1 is taken as default. 

The full renormalisation scale dependence of ([7]) is given by replacing the coupling constant, 
the fixed-order coefficients, the resummation functions and the matching coefficients as follows: 

a s -> a s (fi) , (9) 

B(y) -> B(y, f i) = 2l3 lnx^A(y)+B(y) , 

C(y) -> C(y, f x) = (2p lnx tl ) 2 A(y) + 2lnx fl [2p B(y) + 2p l A(y)}+C(y), (10) 

g 2 (a s L) g 2 (a s L, /x 2 ) = g 2 (a s L) + — (a s L) 2 g[ (a s L) Inx^ , (11) 

G 22 -» G 22 (ji) = G 22 + 20 o G 12 hix ll , 

G 33 -> G33 (a*) = G33 + 4/3 G 23 ln^ . (12) 

In the above, ^ denotes the derivative of g\ with respect to its argument. The LO coefficient 
A and the LL resummation function g±, as well as the matching coefficients Ga+i remain 
independent on [i. 

The arbitrariness in the choice of the logarithm to be resummed can be quantified by varying 
the constant xl ■ This variation implies also the modification of the NLL resummation function 
and of its coefficients 

g 2 {a s L) -> g 2 (a s L) = g 2 {a s L) + — ~ (Lgi(a s L)) lnx L , (13) 
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Gu + 2G 12 \nx L , 
G22 + 3G23I11XL , 
G33 + 4G 34 lnx L . 



(14) 



In Figures [T] and [21 we compare the matched NLLA+NNLO predictions for all event shape 
variables with the fixed order NNLO predictions, and the matched NLLA+NLO with fixed order 
NLO. To allow for a better distinction of the different descriptions, all distributions were weighted 
by the respective shape variables. We use Q = Mz and fix Xn = 1, the strong coupling constant 
is taken as the current world average a s (Mz) = 0.1189 [31]. To quantify the renormalisation 
scale uncertainty, we have varied 1/2 < < 2, resulting in the error band on these figures. 

Several common effects are seen for all shape variables. The most striking observation is 
that the difference between NLLA+NNLO and NNLO is largely restricted to the two-jet region, 
while NLLA+NLO differ in normalisation throughout the full kinematical range. This behaviour 
may serve as a first indication for the numerical smallness of corrections beyond NNLO in the 
three-jet region. 

An immediate consequence of this behaviour concerns the extraction of a s from event shape 
data. Studies at LEP [1-4] yielded substantially different values (by about 10-15%) from NLO 
and NLLA+NLO theory. This discrepancy is an immediate consequence of the varying nor- 
malisations in the two approaches. One can expect that a s obtained using NLLA+NNLO will 
differ from the fixed-order NNLO result [32] only moderately, given the good agreement of both 
descriptions in the three-jet region for fixed a s . 

In the approach to the two-jet region, the NLLA+NLO and NLLA+NNLO predictions agree 
by construction, since the matching suppresses any fixed order terms. Equally, the renormalisa- 
tion scale uncertainty on both these predictions is identical in this region. In the three-jet region, 
NLLA+NNLO agrees with NNLO. The difference between NLLA+NNLO and NLLA+NLO is 
only moderate in the three-jet region, and especially much smaller than the difference between 
the fixed order NNLO and NLO predictions. The renormalisation scale uncertainty in the three- 
jet region is reduced by 20-40% between NLLA+NLO and NLLA+NNLO. 

The parton-level fixed order NNLO and matched NLLA+NLO and NLLA+NNLO predic- 
tions are compared to hadron-level data taken by the ALEPH experiment [1] in Figure [3l The 
description of the hadron-level data improves between parton-level NLLA+NLO and parton- 
level NLLA+NNLO, especially in the three-jet region for most event shapes. The behaviour 
in the two-jet region is described better by the resummed predictions than by the fixed order 
NNLO, although the agreement is far from perfect. This discrepancy was observed already in 
earlier studies based on NLLA+NLO. It can in part be attributed to hadronisation corrections, 
which become large in the approach to the two-jet limit. A very recent study of logarithmic 
corrections beyond NLLA for the thrust distribution [29] also shows that subleading logarithms 
in the two-jet region can account for about half of this discrepancy. 

A precise extraction of a s from event shape data will require the inclusion of hadronisation 
corrections and of quark mass effects (at least to NLO [33]), as done already in the fixed order 
NNLO study [32]. It can be anticipated that inclusion of the matched NLLA+NNLO correc- 
tions results in a further improvement of the extraction of a s from event shape data over results 
obtained previously at NLLA+NLO as well as at NNLO. The principal shortcomings of the 
up-to-now default NLLA+NLO studies were the substantial renormalisation scale uncertainty 
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and the sizable scatter of values of a s obtained from different shape variables. It was observed 
recently, that a fixed-order NNLO extraction [32] reduces the renormalisation scale uncertainty 
by a factor 1.3 compared to NLLA+NLO and eliminates the scatter between different observ- 
ables. It will be very interesting to see the impact of the matched NLLA+NNLO calculation on 
the extraction of a s . We will address this issue in a future study. 

A routine implementing the matching for all event shapes discussed here can be obtained 
upon request from the authors. 
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A Resummation Functions 

We summarize here the expressions for the resummed NLL integrated cross section ([5]) for 
different event shapes. One has 

R(y,Q,ti) = (1 + Cia s ) £ (y) , 

with 

£ (y) = exp {Lgi {a s L) + g 2 (a s L)} . 

Following [22,25], and in order to unify the notation, the resummed part is then expressed 
through auxiliary functions h\ (A) and h 2 (A), with: 

= V s {y)F{R') 

where 

R'{\) = -\[h l {\) + \h! l {\)]. 

The functions h\ (A), h 2 (A), S s (y) and T (R') depend on the event shape observable, as well as 
the parameter A. The QCD constants (3q, (3\ and K are normalised as follows: 

A) = ^ {HO A - 2N F ) , 
Pi = ^ (l7 C\- 5 C A N F -3C F N F 

K = C A («L-£)-*N F . 
I 18 6 / 9 
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A.l Thrust and C-Parameter 

From [19] and [23], one has: 

A = — asL , 

7T 

hl W = -^-[(l-2A)ln(l-2A)-2(l-A)ln(l-A)] , 
h2 (A) = -%^[21n(l-A)-ln(l-2A)]-^ln(l-A) 



( F/ J- fin (1 - 2A) - 2 In (1 - A) + Jin 2 (1 - 2A) - In 2 (1 - A) 

S s (y) = e L " 



2/3 3 V ' ; ' 2 

D L2h 1 (\)+2h 2 (X) 



These yield: 



r(i + 4i?') 



g 1 (a s L) = 2h 1 (^-a s L S j , 

g 2 (a s L) = 2h 2 (^-a s L) -\n[r{l + 4R')]-2 7E R' . 



A. 2 Heavy Jet Mass 

From [19] one has: 

X ft T 

A = — asL , 

7T 

= -^-[(l-2A)ln(l-2A)-2(l-A)ln(l-A)] , 
*»W = -%?[21n(l-A)-ln(l-2A)]-^ln(l-A) 



- fin (1 - 2A) - 2 In (1 - A) + Jin 2 (1 - 2A) - In 2 (1 - A) 



These yield: 



e -2-y E R' 

r(i + 2R'f ' 



gi(a s L) = 2h 1 {^-a s Lj , 

g 2 (a s L) = 2h 2 (^ a s L^j - 2 In [r (1 + 2i?')] - 2 7 £# . 



6 



A. 3 Total Jet Broadening 

From [21,22] one has: 

A 

hi (A) 



9 A) T 

2 — a s L , 

IT 

2C F 



(ln(l-A) + A) , 



h 2 (A) 
S s (y) 



A/3o 
C F K 

+ ^rU ln (1 - A)+ 



(Mi-A) + T ^)-^Mi-A) 



In (1-A) A 



1-A 



1 - A 



_ p Lh 1 {X)+h 2 {\) 



00 dx ( 1 + x 
i x 2 



—R' 



e -2~f E R' 

T(l + 2R') 



4 R ' 2^^, 1 + 2 + R>; -1) 
(l + #) 



r(i + 2 J R / ) 



These yield: 



91 (asL) 

92 (asL) 



h i ( ^ a s L 



TV 



a s L) - In [r(l + 2#)] -2 lE R' 



+2 In 



4 R/ 2 Fi(R', l + R';2 + R>; -1) 
(1 + 1?') 



A. 4 Wide Jet Broadening 

From [21,22] one has: 

A 



9 & r 
2 — a s L , 

7T 



hi (A) 
/» 2 (A) 

S s (y) 
T(R') 



2C F 



(ln(l-A) + A) 



A/? 

Cftf A „ A 
-^( ln(1 - A) + I3A 

2 



3C F 



ln(l - A) 



+^r^ I ^ln 2 (1 - A) + 



L/ii(A)+h 2 (A) 



In (1 - A) A 



1-A 



1-A 



00 /l + x\~ Rl 

1 X 2 



2 e -2 lE R' 

r(i + R'f 



A R ' 2 Fi(R', l + R';2 + R'- -1) 
{l + Ri) 



r(i + i?') 2 



These yield: 



91 {ol s L) 

92 (atsL) 



I' i ( ^7 a s L 

h2 ( — a s L ) - 2 In [r (1 + R')] - 2 lE R' 



7T 



+2 In 



A R ' 2 F 1 {R', 1 + 12'; 2 + 12'; -1) 
(1 + 12') 



A. 5 Two-to-three Jet Transition in the Durham Algorithm 

From [25, 26] one has: 



A 

hi{\) 
h 2 (A) 

S s (y) 



7T 

3C F 



(ln(l-A) + A) 



2/3o 



ln(l-A) 



C F K 



2[5l (1 - A) 



(A + (1-A)ln(l-A)) + 



+ 



C^A + ln(l-A) + l ln2(1 _ A)| 







1-A 



,Lhi{X)+h 2 (X) 



The function T{R') for I3 is known only numerically [25,26]. We interpolate the points using 
a slightly modified version of Newton's divided difference formula implemented in the CERN 
Computer Program Library. These yield: 



91 (atsL) 

92 (asL) 



J 1 ! ( — a s L 



/3b 



h 2 [^a s L) +ln 
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Thrust: y = r = 1 — T and C-parameter: y = C/6 



Gn = 3C F 

G12 = — 2G F 

G22 = ±C F (-169CU + 22N F + 12{C A -4G f )tt 2 ) 

G23 = \C F (-11C A +2N F ) 

G 33 = yig C F [-612C 2 + 180C A JV F + 108G F JV F + (11G A - 27V F ) (-235 C A + 34 N F + 12 (C A - 6 C F ) it 2 ) + 2304 C 2 C(3)] 

G34 = -j^C F (UC A -2N F ) 2 

Heavy jet mass: y = p 

Gn — 3 Cp 

G12 = -2C F 

G22 = ±C F (-imC A +22N F +l2{C A ~2C F )n 2 ) 

G23 = \C F (~llC A +2N F ) 

G33 = j^C F [-612G 2 + 180Ga N f + 108C F N F + (11C A - 2N F ) (-235 C A +34N F + 12 (C A - 3 C F ) tt 2 ) + 576 G 2 C(3)] 



G34 = 


- j^C F (11 C A - 2 N F f 


Total jet broadening: y = B F 


Gn = 
G12 = 
G22 = 
G23 = 
G33 = 

G34 = 


6G F 
-4C* F 

-±C F (35G A - 2N F - 6C A n 2 + 24C F vr 2 + 288C F ln 2 2) 
-§C f (11Ca -2N F ) 

^-G F (-2471G 2 + 7mC A N F + 108G F JV F - 447V 2 + 132G 2 tt 2 - 792G a G f tt 2 - 24G A A r F 7r 2 + 144G F iV F 7r 2 

+864C 2 7r 2 ln2 - 9504G A G F ln 2 2 + 1728G F iV F ln 2 2 - 5184G 2 ln 3 2 + 2376C 2 C (3)) 
-|C F (llC A -2iV F ) 2 


Wide jet broadening: y = Bw 


Gn = 
G12 = 
G22 = 
G23 = 
G33 = 

G34 = 


6G F 
— 4G F 

-§C F (35C A - 2N F - 6C a tt 2 + 288G F ln 2 2) 
-%C F {UC A -2N F ) 

£C F (-2471G A + 760C A N F + 108C F N F - UN 2 + 132G a tt 2 - 2AC A N F ir 2 + 864C 2 7r 2 ln2 - 9504G A G F ln 2 2 

+1728G F JV F ln 2 2 - 5184C 2 ln 3 2 - 2808C 2 C (3)) 
-|C F (11G A -2Af F ) 2 



Two-to-three jet transition in Durham algorithm: y = Y3 



Gn 


= 3C F 


G12 


= -C F 


G22 


= ±C F (-35Ga + 144C F ^2 + 2N F + %C a tt 2 ) 


G23 


= -\C F (UC A -2N F ) 


G33 


= 354 G F (-2471C 2 + 4752C A G F j lr 2 + 2592G 2 T3 + 760C A N F + 108G F JV F - 864G F J c " 2 A f F - 44JV 2 + 132C* 2 n 2 
-24C A N F w 2 ) 


G34 


= -±C F (11C A -2N F ) 2 



Table 1: The logarithmic coefficients Gij for LL and NLL up to the third order in 05. 
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Gu 


Gl2 


G22 


C23 


G33 


C34 


r/C 


4.0 


-2.66667 


-24.9388 


-10.2222 


-285.055 


-45.716 


P 


4.0 


-2.66667 


-13.2415 


-10.2222 


-196.125 


-45.716 


Bf 


8.0 


-5.33333 


-61.8768 


-27.2593 


-824.787 


-156.741 


B w 


8.0 


-5.33333 


-15.0876 


-27.2593 


-472.065 


-156.741 


Y 3 


4.0 


-1.33333 


0.867972 


-3.40741 


-28.1784 


-9.7963 



Table 2: The numerical value of the logarithmic coefficients Gij for LL and NLL up to the third 
order in as- 
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Figure 1: Comparison of the matched NLLA+NNLO and NLLA+NLO with fixed order NNLO 
and NLO predictions for the thrustlike observables r, p and C-parameter. 
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Figure 2: Comparison of the matched NLLA+NNLO and NLLA+NLO with fixed order NNLO 
and NLO predictions for Bt, B\y and Y3. 
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Figure 3: Comparison of the matched NLLA+NNLO and NLLA+NLO with fixed order NNLO 
with the hadron-level data taken by the ALEPH experiment [1] . 
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